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Introduction 


There  exists  a  need  for  innovative  research  to  develop  new  imaging  technologies 
for  sensing  and  quantitating  molecular  signatures  in  vivo  and  in  deep  tissues  The  main 
rationale  for  developing  such  technologies  is  2-fold;  a)  to  allow  the  earlier  detection  of 
disease  based  on  molecular  abnormalities  and  b)  gain  a  better  understanding  of  disease 
mechanisms  to  allow  efficient  targeting  of  specific  signatures  for  early  detection  and  for 
treatment.  The  non-invasive  three-dimensional  mapping  and  quantification  of  cellular  and 
sub-cellular  events  in  living  organisms,  represents  a  highly  novel  and  unexplored  field  in 
imaging  sciences  that  is  expected  to  have  a  major  impact  in  biomedical  research  and  the 
clinical  practice  of  diagnostic  imaging.  The  development  of  novel  optical  tomographic 
techniques  and  molecular  fluorescent  probes  in  the  NIR  can  significantly  enhance 
molecular  and  cell  biology  techniques,  in  the  in-vivo  study  of  cancers  in  deep  tissues. 
Technologies  that  use  light  have  been  developed  in  the  past  to  target  oxy-  and  deoxy¬ 
hemoglobin  or  vascular  contrast  agents.  Under  this  concept  award  we  investigated 
whether  this  technology,  after  necessary  improvements  and  modifications,  would  be 
adequate  to  sense,  reconstruct  and  quantify  fluorescent  reporters  from  within  tissues  in- 
vivo  and  in  a  non-invasive  manner.  Fluorescence-mediated  molecular  tomography  (FMT) 
is  a  sensitive  imaging  modality  that  reconstructs  and  quantifies  the  distribution  of 
fluorescent  reporters  within  tissue  with  high  sensitivity.  The  specific  reporters  targeted 
were  activatable  near-infrared  probes  that  are  non-fluorescing  in  their  native  state  but  are 
activated  in  the  presence  of  specific  enzymes  that  are  over-expressed  in  carcinomas 
For  imaging  purposes,  we  developed  a  highly  sensitive  imager  with  three-dimensional 
imaging  capabilities,  able  to  detect  signals  at  both  intrinsic  and  fluorescent  wavelengths 
We  further  developed  the  necessary  algorithms  for  very  accurate  reconstructions, 
adapted  specifically  to  the  requirements  of  imaging  living  tissue  We  tested  the 
technology  with  phantoms  and  subsequently  applied  it  in-vivo  imaging  of  animal  models. 
The  importance  of  this  imaging  strategy  is  further  amplified  by  considering  that  one  of  the 
major  drawbacks  of  CT  and  MRI  contrast  agents  has  been  the  need  to  deliver  them  at 
pmolar  tissue  concentrations.  Similar  to  nuclear  imaging,  photon  technology  detection 
has  superior  sensitivity  detecting  fluorescent  molecules  at  nano-  to  picomolar 
concentrations.  One  significant  advantage  of  fluorescent  molecules  over  isotopes  is  that 
they  can  be  quenched  and  de-quenched,  enabling  the  design  of  molecular  switches  or 
‘beacons”.  A  further  benefit  is  the  use  of  non-ionizing  radiation  and  the  relatively  low  cost 
instrumentation  required. 


Body 

Imaging  molecular  targets  in  vivo 

The  goal  of  the  research  performed  lies  in  the  development  of  new  sets  of 
strategies  to  interrogate  targets  at  subcellular  levels  (Fig.  1).  These  targets  are  ultimately 
responsible  for  pathogenesis  and  represent  the  key  points  for  targeting  early  detection 
and  subsequent  therapeutic  intervention.  The  impact  of  developing  molecular  imaging 
techniques  is  potentially  enormous.  First,  relatively  crude  parameters  of  tumor  growth  and 
development  (tumor  burden,  anatomic  location  etc.)  could  be  supplemented  by  more 
specific  parameters  (e.g.  growth  factor  receptors,  cell  surface  markers,  transduction 
signals,  cell  cycle  proteins  etc.).  Second,  this  added  information  is  expected  to  go  hand- 
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in-hand  with  the  development  of  novel  targeted  detection  strategies  that  will  allow  us  to 
assess  tumors  at  a  molecular  level,  long  before  phenotypic  changes  occur.  This,  in  turn, 
is  expected  to  have  an  impact  in  the  level  of  specificity  in  detecting  not  only  the  tumor 
presence  but  its  molecular  synthesis  and  further  initiating  appropriate  targeted  treatment 
in  case  of  malignancy.  Third,  molecular  imaging  potentially  allows  one  to  study  the 
genesis  of  diseases  in  the  intact  microenvironment  of  living  systems,  which  is  currently 
very  difficult  to  do  otherwise.  The  importance  of  intact  micro-  and  macroenvironments  on 
tumor  progression  has  recently  been  recognized  Fourth,  molecular  imaging  will  be 
critical  in  testing  novel  drug  delivery  strategies  as  there  exist  significant  barriers  to  the 
delivery  of  complex  and  molecular  therapies^.  Understanding  these  barriers  and  finding 
solutions  to  modulate  them  will  be  of  utmost  importance  in  future  drug  development. 
Finally,  molecular  imaging  allows  one  to  gain  three-dimensional  information  that  may  be 
much  faster  to  obtain  than  is 
currently  possible  with  time 
consuming  and  labor  intensive 
conventional  technique.  The 
key  elements  to  sampling 
molecular  information  are 
summarized  in  Fig.  1.  We 
believe  that  FMT  imaging  with 
activatable  “smart”  imaging 
probes  will  become  a  key 
enabling  technology  towards 

the  above  goals.  F'S  1 :  Comparison  of  in  vitro  and  in  vivo  molecular  imaging. 


“Smart”  imaging  probes 


Photon  detection  has  high  intrinsic  sensitivity  allowing 
the  detection  of  fluorescent  molecules  at  pico-  to  femto- 
molar  concentrations  (see  appendix),  similar  to  isotope 
concentrations  in  nuclear  medicine.  One  significant 
advantage  of  fluorescent  molecules  over  isotopes  is  that 
they  can  be  quenched  and  de-quenched,  enabling  the 
design  of  molecular  switches  or  ‘beacons".  The  underpinning 
mechanism  of  action  of  the  previously  conceived  molecular 
probes  is  illustrated  in  Fig.  2.  The  probes  are  essentially 
non-fluorescent  in  their  native  (quenched)  state  and  become 
highly  fluorescent  after  target  interaction,  resulting  in  signal 
amplification  of  up  to  several  hundred-fold,  depending  on  the 
specific  design.  This  approach  has  four  major  advantages 
over  other  methods  when  single  fluorochromes  are  attached 
to  affinity  molecules:  1)  a  single  enzyme  can  cleave  multiple 
fluorochromes,  thus  resulting  in  one  form  of  signal 
amplification,  2)  reduction  of  background  “noise”  by  several 
orders  of  magnitude  is  possible,  3)  very  specific  enzyme 


Fig.  2:  Example  of  activatable 
enzyme  specific  imaging  probes. 
Note  that  the  fluorochromes  are 
excited  in  the  near  infrared. 


activities  can  potentially  be  interrogated  and  4)  multiple  probes  can  be  arranged  on 


delivery  systems  to  simultaneously  probe  for  a  spectrum  of  enzymes. 
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A  number  of  different  designs  of  optical  smart  agents  have  been  described.  In  this 
work  we  employed  a  probe  designed  to  report  on  specific  enzyme  activities.  The  probes 
typically  consist  of  three  building  blocks:  1)  near  infrared  fluorochromes,  2)  peptide 
substrates  and  3)  a  delivery  vehicle.  The  design  of  the  probes  has  been  described  in 
detail  elsewhere  ^  and  basically  uses  certain  peptide  substrates  (Table  2)  attached  to  a 
delivery  molecule  to  impart  molecular  specificity.  When  a  given  target  protease  cleaves 
the  substrate,  the  fluorochromes  are  released  (see  Fig.2)  and  fluorescence  is  emitted  (de- 
guenching).  Specific  examples  are  tumor  Table 2:  Peptlde  substrates  synthesized  m  our  lab  (The  dots  indicate 

lysosomal  endopeptidases  that  recognize  the  cleavage  site) 

and  cleave  free  lysine  residues  (cathepsin  B) 
cathepsin  D  which  cleaves  two 
phenylalanine  peptides  MMP-2  specific 


probes  or  caspase-3  specific  probes. 


Protease  target 

Peptide  sequence 

Cathepsin  B 

GRR*G 

Cathepsin  D 

GPIC(Et)F*FRLG 

Matrix  metalloproteinase  2 

GPLG*VRG 

Caspase  3 

DEVD'G 

Prostate  specific  antigen 

HSSKLO'G 

Prototype  FMT  imaging  system 

According  to  the  original  statement  of  work  a  prototype  fluorescence-mediated 
molecular  tomography  system  for  3D  quantitative  imaging  in  mice  was  build.  The  system 
components  and  reconstruction  algorithm  were  specifically  chosen  and/or  developed  to  a) 
detect  fluorescence  activation  in  vivo  (i.e.  record  the  smart  probes),  b)  be  quantitative  and 
c)  be  adaptable  to  larger  scales  or  different  geometries.  The  following  is  a  description  of 
the  system  components  while  subsequent  paragraphs  detail  the  experimental  results 
obtained. 


The  FMT  system  consists  of  the  following  components  (Fig.  4  and  Fig.5):  The  light 
source  (a)  is  a  670nm  laser  diode  (B&W  TEK  Newark  DE)  operating  in  continuous  wave 
mode  with  maximum  power  200mW.  This  wavelength  was  selected  for  imaging  the  new 
infrared  fluorochrome  Cy5.5  (Amersham  Pharmacia  biotech  Piscataway  NJ),  which  has 
an  absorption  peak  at  675nm  and  an  emission  peak  at  694nm.  Laser  diodes  at  different 
wavelengths  could  be  used  to  excite  different  fluorochromes.  The  light  from  the  laser 
diode  is  directed  to  a  50/50  fused  silica  coupler  (b)  (OZ  Optics  Ontarion  Canada)  that 
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splits  the  signal  to  the  reference  branch  (c)  and  to  the  optical  switch  (d).  The  optical 
switch  (Dicon  FiberOptics  CA)  can  multiplex  one  input  to  any  of  the  32  outputs.  In  the 
current  implementation  light  is  switched  sequentially  to  only  24  outputs  and  directed  onto 
the  optical  bore  (f)  using  24  200/230  Dm  multimode  optical  fibers  (e)  (Fiber  Instrument 
Sales  Oriskany  NY).  Three-millimeter  diameter  optical  fiber  bundles  (Dollan-Jenner 
Industries  MA)  are  used  for  light  detection  fg).  The  detection  fibers  are  arranged  in  a  6x6 
array  onto  a  flat  plate  (h)  made  of  Delrin™.  The  fiber  array  is  imaged  with  a  Vers  Array 
51 2B  back-illuminated  CCD  camera  (Roper  Scientific  NJ)  with  ~80%  quantum  efficiency 
in  the  700nm  spectral  range. 

The  imaging  lens  is  a  60mm 
F-mount  Micro-Nikkor  f/2.8D 
lens  (Nikon  JP).  For 
fluorescence  detection  a 
three-cavity  band-pass 

interference  filter  is  used  with 
center  frequency  at  705  nm 
and  a  full-width  at  half¬ 
maximum  response  of  +5nm 
(Andover  Corporation  NH). 

Fig.  5:  System  schematic  (see  text  for  details) 

The  optical  bore  (f)  is 

the  central  piece  of  the  device.  It  is  made  out  of  Delrin  and  has  an  inner  diameter  of 
2.5cm  and  an  outer  diameter  of  5cm.  The  optical  fibers  are  inserted  through  the  Delrin 
and  come  face  to  face  with  the  inner  wall.  Three  rings  of  12  detector  fibers  are  interlaced 
with  2  rings  of  12  source  fibers.  The  distance  between  rings  is  0.5cm.  The  selected  bore 
diameter  can  comfortably  accommodate  most  mouse  strains.  The  tube  is  water-sealed 
and  is  filled  with  a  diffuse  matching  fluid  to  minimize  photon  wave  mismatches  within  the 
cylinder.  The  matching  fluid  consists  of  water,  an  absorbing  dye  and  Intralipid  or  Ti02 
particles  at  concentrations  optimized  to  simulate  tissue  absorption  and  scattering. 
Phantoms  and  animals  are  inserted  in  the  optical  chamber.  The  selected  diameter  of  2.5 
cm  can  accommodate  all  mouse  models  proposed. 


Analytical  instrument  evaluation  is  included  in  the  appendix  (report  #1). 


Novel  reconstruction  algorithm 

In  addition  to  developing  a  system  we  also  developed  an  appropriate  self- 
calibrated  algorithm  for  FMT  imaging  based  on  the  general  mainframe  of  perturbation 
algorithms  previously  developed  for  DOT  The  specific  algorithm  developed 
capitalizes  on  the  instrument  ability  to  obtain  intrinsic  and  fluorescence  measurements 
and  to  construct  differential  measurements  not  using  pre-/post-  probe  administration 
measurements  but  using  absorption/fluorescence  differential  measurements.  Specifically, 
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the  intrinsic  photon  field  divides  the  fluorescence  measurement  to  normalize  it  and  to 
cancel  out  instrument  gain  terms  as  well.  In  a  manuscript  recently  published  ^  (see  also 
appendix)  we  show  the  construction  of  the  normalized  field  and  of  appropriate  Greens 
functions  and  demonstrate  that  one  could  use  only  post-administration  measurements  to 
perform  differential  imaging.  An  advantage  of  this  approach  is  that  no  instrument 
calibration  is  required.  Furthermore,  the  algorithm  developed  is  suited  for  imaging 
heterogeneous  media  because  it  has  been  shown  that  perturbative  diffuse  optical 
tomography  methods  perform  very  well  when  differential  measurements  are  obtained 
Although  the  method  is  based  on  a  linear  analytical  solution,  the  basic  principle  of  using 
the  intrinsic  field  to  normalize  fluorescence  measurements  can  be  used  with  higher  order 
analytical  or  numerical  solutions  of  the  diffusion  equation. 

Analytical  algorithm  description  and  evaluation  is  included  in  the  appendix  (report 

#2). 


Phantom  measurements:  calibration,  quantitation,  validation,  sensitivity 


Absorption  Imaging.  Phantom  experiments  were  initially  performed  to  verify  the  three- 
dimensional  position  and  accuracy  of  measurements  in  both  absorption  and  fluorescence 
mode.  The  experimental  set-up  is  illustrated  in  Fig  6.  Briefly,  a  phantom  containing  3 
capillary  tubes  (1  mm  internal 
diameter),  was  constructed  using  a 
triangular  geometry  and  inserted 
into  the  optical  chamber  containing 
a  turbid  medium  (Intralipid).  The 
capillaries  were  separated  8  and 
1 1  mm  from  each  other,  and  were 
coated  with  a  black  dye  to 
maximize  absorption.  Fig.  6  depicts 
the  reconstructed  image  at  a  plane 
perpendicular  to  the  longitudinal 
axis  of  the  imaging  chamber,  at 
about  the  middle  of  the  three- 
dimensional  volume.  The 
reconstruction  mesh  used  was  0.8  x  0.8  x  2  mm^.  The  reconstruction  used  24  x  36 
measurements.  Our  data  show  that  very  good  positional  accuracy  can  be  obtained  using 
the  standard  perturbation  based  on  the  Rytov  approximation.  For  this  measurement  the 
Greens  functions  are  derived  using  a  finite-difference  solution  of  the  diffusion  equation 
assuming  a  homogeneous  cylindrical  medium  with  the  Intralipid  optical  properties  at  675 
nm  and  bounded  by  a  high  absorption  cylindrical  boundary  (pa=5cm'’). 


A.  Phantom  set-up 


B.  Reconstructed  image 


Fig,  6:  Expeiimetal  set-up  and  actual  reconstruction  (right).  Note  the 
excellent  spatial  location  of  the  capillary  tubes 
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Sensitivity.  Subsequent 
phantom  experiments  were 
performed  to  investigate  the 
sensitivity  of  detection  using 
non-quenched  near  infrared 
Cy5.5  dye.  We  inserted  a  3 
mm  diameter,  2  cm  long  tube 
containing  200nM  Cy5.5  in  the 
optical  bore  (Fig.  7).  The 
background  intralipid  solution 
had  optical  properties: 

[Xa=0.025  cm'^  and  ps-S  cm'\ 

The  total  fluorescent  photon 
counts  detected  from  two 
opposite  sides  after  noise  subtraction  are  plotted  as  a  function  of  light  power.  The  light 
source  was  adjacent  to  position  A  and  the  integration  time  was  20  seconds.  The  data 
show  that  there  is  significant  fluorescent  signal  detected  at  all  positions.  Increasing  the 
laser  power  can  significantly  reduce  acquisition  times.  Based  on  those  measurements 
and  appropriate  solutions  of  the  diffusion  equation  we  can  predict  the  equivalent  photon 
counts  expected  in  different  geometries  and  tissues.  For  example  within  a  slightly 
compressed  breast  of  5  cm  width  and  optical  properties  pa=0.05cm'^  and  ps’=10cm'^  the 
expected  photon  counts  from  a  similar  3  mm  volume  and  expected  dye  concentrations 
the  maximum  counts  expected  would  be  more  than  2,500  above  noise  level  using  a  100 
mW  laser  power.  Other  measurements  in  the  laboratory  with  fluorometers  using  higher 
light  powers  suggest  that  the  signal  detected  should  scale  linearly  with  light  power  for 
much  higher  powers  used. 


Fig.  7:  Sensitivity  of  read-out  of  fluorochrome  phantom  in  Intralipid. 


[Cy5.5] 
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Quantitating  fluorochrome  concentration  by  imaging. 

To  determine  the  accuracy  in  estimating 
fluorochrome  concentrations  we  also  imaged  a  similar 
phantom  as  described  above  but  containing  different 
concentrations.  Again,  the  optical  bore  was  filled  with 
0.5%  Intralipid  (Fresenius  Kabi  Clayton  L.P.,  NC)  and 
6pM  IC-Green  (Akorn  Inc.  IL)  yielding  a  liquid  with  ^i[=^ 
cm'^  and  =0.05  cm'^  at  675nm.  Two  glass  tubes  (3.5 

mm  inner  diameter)  were  inserted  with  centers  12  mm 
apart.  The  tube  on  the  left  contained  500  nM  Cy5.5  and 
the  tube  contained  750  nM  Cy5.5.  The  reconstructed 
plane  was  interpolated  to  a  60  x  60  grid  (h’=  0.04  x  0.04 
cm^).  Two  objects  were  reconstructed  and  the  accuracy 
of  quantitating  the  true  Cy5.5  concentration  was  ±  10%. 

The  result  is  shown  in  Fig,  8  and  analytically  in  the  appendix  (report  #2).  The 
reconstruction  accuracy  can  be  further  improved  with  larger  data  sets 


Fig.  8:  Reconstructions  of  2  phantoms 
containing  different  concentrations  of  Cy 
5.5 
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Imaging  of  fluorescence  activation. 

In  another  experiment,  quantitative,  spatially  localized  information  on  fluorescence 
activation  was  obtained  as  a  function  of  time.  As  shown  in  Fig.  9,  a  3  mm  tube  was 
immersed  into  Intralipid,  similarly  as  above.  The  tube  contained  1.5  pM  of  a  trypsin 
sensitive  enzyme  probe,  which  was  activated  by  trypsin  (100  pM)  at  time  0.  Each  single 
frame  was  acquired  for  1  minute  (12  sources  at  5  sec  acquisition  each).  The  series  of 
temporal  snapshots  on  the  right  of  Fig.  9  illustrate  the  reconstructed  frames  obtained  at 
different  time  points.  For  example,  at  20  minutes  after  trypsin  addition,  ~20%  of  the  probe 
had  been  activated.  At  later  time  points,  e.g.  at  200  minutes  after  the  enzyme  was  added, 
~75%  of  the  probe  had  been 
activated.  These  experimental 
data  show  that  the  developed 
FMT  system  can  monitor  with 
significant  sensitivity  and 
accuracy  fluorescence 

activation  transient  phenomena 
as  corroborated  by  parallel  in 
vitro  enzymologic 

measurements.  It  is  expected, 
that  the  proposed 
improvements  to  the  FMT 
system  will  1)  further  improve 
the  sensitivity  by  improving  the 
light  strength,  2)  improve 
quantification  accuracy  but 
improving  laser  intensity 
stability  and  3)  improve 
acquisition  speeds  by 

appropriate  automation  and  the  Fig.  9:  Temporal  imaging  of  enzyme  activity  in  diffuse  media  (see  text) 

proposed  power  increase. 


In-vivo  imaging 

We  tested  the  novel  imager  with  animal  imaging  as  well.  To  perform  validation  we 
combined  FMT  measurements  with  anatomical  and  gadolinium  enhanced  MRI.  For 
coregistration  purposes  the  animals  were  placed  in  an  appropriately  constructed 
translucent  insert  (shown  in  Fig.  10)  that  was  transferable  between  the  optical  and  MR 
scanners  to  maintain  positional  certainty.  We  used  this  combined  MRI  and  FMT  scan  to 
obtain  maps  of  cathepsin  B  protease  activity  in  human  tumors  implanted  into  nude  mice 
(Fig.  11).  The  tumors  were  cathepsin  B  rich  HT1080  fibrosarcoma,  which  had  been 
implanted  into  the  mammary  fat  pad  of  mice  7-10  days  prior  to  the  experiment.  The 
animals  received  an  IV  injection  of  a  cathepsin  B  sensitive  imaging  probe  [Weissleder, 
1999  #305]  at  24  hours  prior  to  the  imaging  experiments.  The  animals  were 
anaesthetized  with  an  intraperitoneal  injection  of  90  mg/kg  ketamine  and  9  mg/kg  xylazine 
and  were  imaged  by  FMT  (manual  image  acquisition  was  approximately  10  minutes). 
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Subsequently,  the 
animals  were  removed 
from  the  imaging 
chamber.  Fiducials  with 
water  were  attached  to 
predetermined 
positions  on  the 
periphery  of  the  insert. 

The  insert  was 
subsequently  placed  in 
the  MR  coil  and  a  set 
of  axial  T2-weighted 
imaged  were  obtained. 

The  role  of  the  fiducials  was  to  identify  on  the  MR  images  the  position  of  selected  source 
and  detector  fibers  for  later  co-registration  of  the  images.  Fig.  9  shows  three 
corresponding  slices  obtained  at  tumor,  heart  and  kidney  levels.  There  was  excellent 
congruence  of  optical  and  MR  contrast  from  the  tumor  slice,  known  to  have  high 
cathepsin  B  activity  as  corroborated  by  immunohistochemistry  and  Western  blotting.  The 
cardiac/pulmonary  slice  (not  shown  here)  demonstrated  cathepsin  B  absence  in  those 
organs  and  the  lowest  slice  (not  shown  here)  shows  renal  presence  of  the  enzyme  as 
reported  in  the  literature. 


Insert  for  FMT  imaging  Same  inserts  for  MRI  coil  Object  in  MRI  coil 


Fig.  10:  Various  mouse  holders  (inserts)  for  the  FMT  imaging  system  (left).  The  same  inserts 
also  fit  snuggly  into  home  built-MR  coils  for  combined  MR/FMT  imaging. 


MRI  Image  FMT  Image  Merged  Image 


0  0.2  0.4  0.6  0.8  1.0  1.2 


Fig.  11:  Combined  MR/FMT  measurement  at  three  body  levels  of  a  tumor  bearing  mouse.  FMT  measure  the 
cathepsin  B  enzyme  activity 


To  test  whether  the  system  and  technique  could  investigate  smaller  tumors, 
implanted  in  deeper  sites  we  selected  a  brain  implantation  site  and  imaged  millimeter 
sized  tumors  at  different  depths.  A  report  on  that  can  be  seen  in  the  appendix  (report  #3). 
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Key  Research  Accomplishments 

The  key  research  accomplishments  associated  with  the  statement  of  work  are  as  follows: 

1)  Construction  of  a  fully  operational  FMT  scanner  prototype  for  animal  imaging.  The 
specifics  of  the  instrument  are  included  in  the  appendix  (report  #1). 

2)  Development  of  appropriate  algorithms  for  in-vivo  imaging  that  produce  very 
accurate  quantified  images  of  the  activation  of  smart  fluorescent  probes  and  molecular 
beacons. 

3)  Testing  of  the  above  instrument  with  phantom  measurements  to  demonstrate 
femto-molar  sensitivity  to  detect  fluorescent  probes  [see  appendix  report  2]  and  the 
capacity  to  perform  three-dimensional  quantified  imaging  of  fluorochromes  in  dense 
media  [see  appendix  report  1]. 

4)  In-vivo  demonstration  of  non-invasively  imaging  cathepsin  B  activity  in  different  sized 
tumors  implanted  at  different  depths,  correlated  with  MRI,  immuno-histochemistry  and 
Western  blotting. 

5)  Proof  of  feasibility  that  Fluorescence-mediated  molecular  tomography  can  be  used  as  a 
novel  workframe  for  examine  molecular  events  in-vivo  and  in  deep  tissues. 


Reportable  outcomes 


1)  Paper  in  preparation  regarding  the  instrument  specifics. 

2)  Paper  in  preparation  regarding  in-vivo  imaging. 

3)  Paper  published  (Ntziachristos  et.al.  Opt.  Lett.  26(12):  893-895  (2001))  on 
algorithmic  development. 

4)  Presentation  on  FMT  in  High  resolution  in  small  animals  conference  2001 
(HiRes2001)  held  in  Rockville  Maryland,  September  2001. 

5)  Invited  presentation  in  Contrast  Media  Research  2001  (CMR  2001)  held  in  Capri, 
Italy,  October  2001. 

6)  R21/R33  grant  submitted  to  NIH  for  continuation  of  instrument  development. 


Conclusions 

We  have  developed  an  instrument  and  an  imaging  method  and  applied 
successfully  to  imaging  cancer-associated  expression  of  cathepsin  B  probes.  We  found 
that  the  technique  offeree  suberb  sensitivity,  quantification  and  localization  accuracy. 
Clearly  higher  source  detector  systems  could  further  improve  the  sensitivity,  accuracy  and 
resolution  which  in  the  current  implantation  is  of  the  order  of  2  mm.  The  technique  could 
serve  as  the  work  platform  for  testing  multiple  other  cancer  associated  enzymes  or  for 
gene  expression  profiling.  Implementation  of  spectral  capabilities  could  further  enhance 
this  approach  since  multiple  targets  could  be  investigated  simultaneously  by  separating 
their  contribution  using  fluorochromes  emiting  at  different  wavelengths.  The  former  is  of 
particular  clinical  relevance  as  “expression  patterns”  rather  than  single  molecules  may 
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form  the  basis  for  more  sophisticated  future  imaging  diagnostic  techniques.  Most 
importantly  we  have  found  that  this  technology  could  be  transferable  to  human 
applications.  NIR  light  has  been  shown  to  penetrate  tissues  several  centimeters  and 
detection  is  thus  not  necessarily  limited  to  small  volumes.  Indeed,  recent  measurements 
and  modeling  [Ntziachristos,  2001  #30]  show  that  near  infrared  fluorescent  signals  from 
tumor-like  structures  can  propagate  for  more  than  15  centimeters  in  breast  or  lung  tissue 
and  more  than  5  cm  in  the  adult  brain.  Again,  these  estimates  may  well  be  surpassed  by 
newer  detection  technologies  and/or  more  efficient  beacons.  While  we  chose  to  construct 
a  cylindrical  imaging  chamber  for  mice,  human  imaging  systems  could  well  have  different 
geometries  for  more  efficient  light  collection.  Planar  imaging  geometries,  handheld 
scanners  or  even  endoscopic  sensors  are  all  technologically  feasible  providing  optimized 
approaches  for  imaging  of  specific  organs  or  diseases.  An  important  aspect  to  potential 
clinical  tomographic  imaging  is  the  fact  that  FMT  is  inherently  quantitative  and  can  be 
acquired  over  time  for  monitoring  purposes.  Additional  advantages  include  the  fact  that  no 
radiation  is  required,  that  beacons  and  fluorochromes  are  usually  stable  and  do  not  decay 
like  isotopes  and  that  the  technology  is  relatively  inexpensive  compared  to  other 
tomographic  imaging  systems.  Finally,  FMT  imaging  should  be  useful  for  monitoring  drug 
therapies  at  the  molecular  levels  [Bremer,  2001  #21]  over  time  and  in  tailoring  treatments 
in  individual  patients.  FMT  imaging  may  also  be  a  useful  adjunct  to  x-ray  computed 
tomography  and  magnetic  resonance  imaging.  It  is  thus  hoped  that  FMT  imaging  will  find 
applications  in  unraveling  complex  biological  pathways  in  murine  models  and  In  applying 
this  knowledge  to  current  problem  areas  in  clinical  diagnosis  and  drug  testing. 
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CCD-based  scanner  for  three-dimensional  fluorescence-mediated  diffuse 
optical  tomography  of  small  animals. 
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Massachusetts  General  Hospital  &  Harvard  Medical  School. 
Charlestown  MA  02129 


ABSTRACT 

We  present  a  novel  tomographer  for  three-dimensional  reconstructions  of  extrinsic 
fluorochromes  distributed  in  small  animals.  The  tomographer  has  been  developed  for  studying  novel 
classes  of  fluorescent  agents  useful  in  diagnostics  and  drug  development.  Photon  detection  is  based 
on  CCD  technology  and  allows  the  implementation  of  a  large  parallel  array  of  optical  channels  with 
high  sensitivity.  Using  this  instrument  we  studied  the  response  and  detection  limits  of  near-infrared 
fluorochromes  in  diffuse  media  as  a  function  of  light  intensity  and  for  a  wide  range  of  biologically 
relevant  concentrations.  We  further  examined  the  resolution  of  the  scanner  and  the  reconstruction 
linearity  achieved  in  the  fluorescence  mode.  We  demonstrate  that  the  instrument  attains  better  than 
3mm  resolution,  is  linear  within  more  than  two  orders  of  magnitude  of  fluorochrome  concentration 
and  can  detect  fluorescent  objects  at  sub-nanomolar  fluorochrome  concentrations. 


INTRODUCTION 

Diffuse  Optical  Tomography  (DOT)  is  a  medical  imaging  techrtique  that  collects  light 
propagated  through  diffuse  media,  such  as  tissue,  at  multiple  projections  and  reconstructs  the 
internal  distribution  of  optical  properties  '.  Interestingly,  the  teclmique  propagated  directly  from  the 
laboratory  bench  to  human  applications,  such  as  imaging  of  breast  cancer  ^  3,  brain  function  ^  and 
muscle  3  because  of  the  non-invasiveness  and  safety  of  light  radiation  used.  It  would  be  of  great 
interest  however  to  develop  instruments  for  tomographic  animal  investigations,  particularly  for  in- 
vivo  studies  of  novel  contrast  agents  and  fluorescent  probes  which  are  not  yet  approved  for  human 
use. 


In  this  work  we  describe  a  tomographic  scanner  operating  in  fluorescence  mode,  developed 
specifically  for  imaging  of  small  rodents.  The  scanner  offers  highly  sensitive  detection  of 
fluorochromes  and  high  capacity  for  parallel  detection  of  a  large  number  o,f  optical  fibers,  in  order  to 
facilitate  true  three-dimensional  reconstmctions.  We  have  recently  reported  on  the  capacity  of  the 
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scanner  to  perform  quantitative  three-dimensional  reconstructions  of  fluorochromes  using  a  newly 
developed  fluorescence  reconstmction  algorithm  Here  we  focus  on  studying  key  important 
experimental  parameters  on  system  performance  and  fluorochrome  reconstmctions,  i.e.  the  detection 
sensitivity  and  the  linearity  of  fluorescent  signals  as  a  function  of  fluorochrome  concentration  and 
incident  light  intensity.  Additionally,  we  examine  the  linearity  achieved  in  the  reconstructed 
fluorescence  images  as  a  function  of  fluorochrome  concentration.  To  the  best  of  our  knowledge  this 
is  the  first  systematic  experimental  evaluation  of  tomographic  performance  of  fluorochromes 
embedded  in  diffuse  media.  To  compete  the  study  we  also  report  on  the  signal  to  noise  and  the 
resolution  achieved. 

The  remainder  of  this  paper  is  divided  in  three  sections.  The  first  section  describes  the 
instrument  developed  and  the  reconstruction  tools  employed.  The  second  section  presents 
experimental  measurements  from  fluorochromes  embedded  into  diffuse  media.  Finally,  the  third 
section  discusses  the  results  and  outlines  the  potential  applications  of  this  technology. 


METHODS 

Fluorescence  animal  tomoprabher 

The  system  developed  is  shown  in  figure  1.  The  light  source  (a)  is  a  670nm  laser  diode 
(B&W  TEK  Newark  DE)  operating  in  continuous  wave  mode  with  maximum  power  200mW.  This 
wavelength  was  selected  for  imaging  the  new  infrared  fluorochrome  Cy5.5  (Amersham  Pharmacia 
biotech  Piscataway  NJ),  which  has  an  absorption  peak  at  675nm  and  an  emission  peak  at  694nm. 
Laser  diodes  at  different  wavelengths  could  be  used  to  excite  different  fluorochromes.  The  light  from 
the  laser  diode  is  directed  to  a  50/50  fused  silica  coupler  (b)  (OZ  Optics  Ontarion  Canada)  that  splits 
the  signal  to  the  reference  branch  (c)  and  to  the  optical  switch  (d).  The  optical  switch  (Dicon 
FiberOptics  CA)  can  multiplex  one  input  to  any  of  the  32  outputs.  In  the  current  implementation 
light  is  switched  sequentially  to  only  24  outputs  and  directed  onto  the  optical  bore  (f)  using  24 
200/230  /tm  multimode  optical  fibers  (e)  (Fiber  Instrument  Sales  Oriskany  NY).  Three-miUimeter 
diameter  optical  fiber  bundles  (DoUan-Jenner  Industries  MA)  are  used  for  light  detection  (g).  The 
detection  fibers  are  arranged  in  a  6x6  array  onto  a  flat  plate  (h)  made  of  Dekin™.  The  fiber  array  is 
imaged  with  a  Vers  Array  512B  back-illuminated  CCD  camera  (Roper  Scientific  NJ)  with  ~80% 
quantum  efficiency  in  the  700nm  spectral  range.  The  imaging  lens  is  a  60mm  F-moimt  Micro-Nikkor 
f/2.8D  lens  (Nikon  JP).  For  fluorescence  detection  a  three-cavity  band-pass  interference  filter  is  used 
with  center  frequency  at  705  nm  and  a  full-width  at  half-maximum  response  of  ±5nm  (Andover 
Corporation  NH). 

The  optical  bore  (f)  is  the  central  piece  of  the  device.  It  is  made  out  of  Deltin  and  has  an 
inner  diameter  of  2.5cm  and  an  outer  diameter  of  5cm.  The  bore  is  shown  in  more  detail  in  figure  2. 
The  optical  fibers  are  inserted  through  the  Delrin  and  come  face  to  face  with  the  inner  wall.  Three 
rings  of  12  detector  fibers  are  interlaced  with  2  rings  of  12  source  fibers.  The  distance  between  rings 
is  0.5cm.  The  selected  bore  diameter  can  comfortably  accommodate  most  mouse  strains.  The  tube  is 
water-sealed  and  is  filled  with  a  diffuse  matching  fluid  to  minimize  photon  wave  mismatches  within 
the  cylinder.  The  matching  fluid  consists  of  water,  an  absorbing  dye  and  IntraHpid  or  Ti02  particles 
at  concentrations  optimized  to  simulate  tissue  absorption  and  scattering. 

Fluorescence  Diffuse  Optical  Tomoprap/n 

In  this  work  we  employed  a  normalized  reconstruction  algorithm  reported  earlier  ^  that 
operates  on  relative  photon  signals,  thus  facilitating  minimum  instrument  calibration  (in  contrast  to 
operating  on  absolute  photon  measurements).  The  algorithm,  briefly  described  here  for  presentation 
completeness,  uses  measurements  that  can  be  obtained  after  dye  administration  and  resolves  absolute 
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fluorochrome  concentrations;  therefore  it  is  ideally  suited  for  in-vivo  measurements.  For  every 
source-detector  pair,  the  normalized  fluorescence  measurement  can  be  constructed  as: 


_1_  QE^' 


(1) 


where  U jj  (r^ ,  )  is  the  fluorescence  field  detected  at  position  due  to  a  source  at  position  and 
^inc (K incident  field,  i.e.  the  photon  field  at  the  excitation  wavelength  for  the  same 
source  detector  pair.  QE^\  QE^^  are  the  detector  quantum  efficiencies  at  wavelength  ^1^X2 
respectively  and  0yis  the  attenuation  of  the  filter  that  is  used  to  collect  the  fluorescent  field.  The 
terms  0^  and  QE^^IQE^'  can  be  determined  experimentally  and  are  virtually  independent  of 
position  and  time.  Therefore  the  term  U"^  (7^,  Pj)  is  independent  of  individual  source-detector  pair 

gain  factors.  It  can  be  shown  (see  ref#'^  and  references  therein)  that  U"^ can  be  written  as  a 
function  of  the  medium’s  fluorochrome  distribution,  i.e.: 


U"\7„7,)  = 


1 


U,{7„7,X') 


(2) 


t»r(F)  D 


where  is  the  diffusion  coefficient  at  Xi,  U  is  the  speed  of  light  into  the  medium  and  n{7)  is  the 
position  dependent  product  of  the  fluorochrome  absorption  coefficient  (?)  and  fluorescence 
quantum  yield  / .  The  parameters  are  the  wave  propagation  vectors  ^  at  wavelengths  X|,  I2 

respectively  and  r(?)  is  the  fluorescence  hfetime  of  the  fluorochrome  at  position  7  .  Eq.2  is  written 
for  the  general  case  where  the  source  term  is  light  of  modulated  intensity  a .  In  the  constant  wave 
case  CO  =0  and  the  fluorescence  hfetime  r(?)  does  not  contribute  to  the  detected  photon  field.  The 

function  U q{7^,7 )  describes  the  estabhshed  photon  field  into  the  diffuse  medium  at  position 

?  due  to  a  source  at  position  7^  and  G(7j—7,k^^)  is  the  Green’s  function  solution  to  the 
diffusion  equation  and  describes  the  propagation  of  the  emission  photon  wave  from  the 
fluorochrome  at  position  7  to  the  detector  at  7^ .  Eq.2  treats  fluorochromes  as  two-level  quantum 

systems  and  ignores  fluorescence  saturation  effects  and  high-order  scattering  interactions.  These 
approximations  are  generally  met  at  low  fluorochrome  concentrations.  We  also  demonstrate 
experimentally  in  the  results  section  that  saturation  effects  are  not  observed  over  more  than  two 
orders  of  magnitude  of  biologically  relevant  fluorochrome  concentrations. 

For  reconstruction  purposes  Eq.  3  is  discretized  into  a  number  of  volume  elements  (voxels), 
which  yields  a  set  of  coupled  linear  equations.  The  terms  Uo  and  G  are  calculated  numerically  or 
analytically  for  the  appropriate  boundary  conditions  and  background  optical  properties  at  A,i  and  Xz 
respectively.  In  this  work  the  functions  Uo  and  G  were  calcxilated  for  the  cylindrical  boundary  using 
a  three-dimensional  finite-difference  solution  of  the  homogenous  diffusion  equation  for  the 
background  medium  with  discretization  step  h=0.5mm.  The  zero-boundary  condition  *  was  assumed 
at  1///’ outside  the  actual  boundary.  The  resulting  system  of  equations  is  then  inverted  for  the 
unknown  quantity  tz(?)  using  the  Algebraic  Reconstruction  Technique  with  relaxation  parameter 
p=Q.\.  Finally  fluorochrome  concentration  [?^(f)]  is  calculated  as  [?’(f )]  =  ^ ,  where  f  is 

the  fluorochrome’s  extinction  coefficient. 
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RESULTS 

This  section  describes  experimental  measurements  that  outline  the  instrument  performance 
and  the  reconstmction  capacity  pertaining  to  quantification  accuracy,  linearity,  sensitivity  and 
resolution. 

Simal  to  noise  ratio 

To  characteri2e  the  signal  to  noise  ratio  (SNR)  of  the  instrument  we  performed 
measurements  of  intrinsic  light  (without  the  fluorescence  filter)  employing  one  source  and  all  36 
detectors.  Signal  intensity  was  adjusted  so  that  maximum  signal  was  obtained  without  saturating  any 
of  the  CCD  pixels.  Then  200  measurements  on  a  diffuse  medium  (/4,=0.02cm-'  //j  — 10cm-')  were 
sequentially  acquired  over  a  period  of  5  minutes.  This  time  is  representative  of  typical  acquisition 
times  at  full  three-dimensional  tomographic  mode.  The  integration  time  per  measurement  was  0.1 
sec.  The  signal  to  noise  ratio  in  decibels  (dB)  was  calculated  as  iSW/?  =  20  •  log]o(5  /<j)  where  S 

was  the  mean  intensity  of  the  detected  signal  T  after  background  subtraction  and  cr  was  its  standard 
deviation.  The  results  for  all  36  detectors  employed  are  shown  in  figure  3.  The  average  SNR  obtained 
was  ~50dB. 

The  main  source  of  noise,  besides  the  expected  Poisson-distributed  photonic  noise  and 
CCD  read  noise,  was  the  intensity  fluctuation  of  the  laser  source.  The  signal  to  noise  ratio  can  be 
improved  using  the  reference  fiber  measurements  to  monitor  laser  fluctuations  as  a  function  of  time 
and  accordingly  correct  the  measurements.  Figure  3  depicts  the  signal  to  noise  ratio  achieved  for  all 
36  detectors  after  data  correction.  Correction  consists  of  dividing  aU  detector  measurements  at  each 
time  point  with  the  reference  measurement  obtained  at  the  same  time  point  and  multiplying  the 
result  with  the  maximum  reference  value  observed  throughout  the  200  measurements.  The  SNR 
improvement  observed  after  correction  is  approximately  7dB,  since  long-  term  drift  and  other 
fluctuations  can  be  effectively  accounted  for. 


Resolution 

The  resolution  of  tomographic  techniques  in  the  diffuse  regime  is  a  function  of  the  number 
of  sources  and  detectors  employed,  their  geometrical  arrangement,  the  overall  dimension  and  optical 
properties  of  the  medium  of  investigation  and  the  signal  to  noise  ratio.  Fiuthermore,  the  resolution 
achieved  is  not  uniform  but  it  decreases  as  a  function  of  depth.  Here,  we  have  investigated  the 
capacity  of  the  scanner  to  resolve  two  fluorophores  at  the  middle  of  the  optical  bore,  which 
represents  a  worst  resolution  case.  We  inserted  two  Pasteur  glass  tubes  at  the  middle  of  the  optical 
bore  as  shown  in  figure  4a.  The  tubes  had  an  outer  diameter  of  ~1.5mm,  an  inner  diameter  of 
~1.2mm  inner  and  were  filled  with  800nM  of  Cy5.5.  The  optical  bote  was  filled  with  Intralipid  and 
ink  and  attained  a  background  absorption  coefficient  //a=0.05cm-’  and  reduced  scattering  coefficient 
Ps’=  8cm-'.  The  two  tubes  were  placed  at  different  distances  to  each  other  and  imaging  was 
performed  to  identify  the  resolution  limit. 

Figure  4b  and  c  depict  the  reconstmctions  achieved  when  the  tubes  were  placed  with  a  clear 
gap  of  3mm  and  2.8  mm  apart.  In  both  cases  the  two  objects  are  resolved  but  at  the  shorter  distance 
symmetry  is  lost.  When  the  objects  were  placed  even  closer  to  each  other,  the  asymmetry  increased 
until  only  a  single  object  was  resolved.  The  appearance  of  asymmetry  is  mainly  due  to  noise,  which 
affects  the  ill-posed  inversion  problem.  It  should  be  noted  that  the  instrument  constructed  is  not 
built  to  favor  resolution  since  only  12  measurements  per  ring  are  obtained.  Increasing  the  number  of 
sources  and  detectors  along  the  periphery  can  further  improve  resolution,  due  to  the  higher  spatial 
sampling  but  also  due  to  the  inherent  improvement  in  signal  to  noise  ratio  Nevertheless,  under  the 
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current  implementation,  the  resolution  of  the  imager  build  is  better  than  3mm,  especially  since  the 
phantom  examined  is  a  worse  resolution  case. 

'Fluorescence  as  a  function  of  concentration  and  lifht  intensip. 

A  key  parameter  of  DOT  in  fluorescence  mode  is  the  emission  signal  response  as  a  function 
of  light  intensity  and  fluorochrome  concentration.  Deviation  from  a  linear  behavior  indicates  photo- 
bleaching  or  quenching  phenomena  that  could  result  in  quantification  errors  and  data 
misinterpretation.  In  order  to  investigate  fluorescence  responses,  we  inserted  a  fused  quartz  tube 
(Wilmad  Glass  NJ)  into  the  optical  bore  as  shown  in  figure  5.  The  tube  was  4mm  in  inner  diameter 
and  8mm  in  height  and  was  sequentially  filled  with  different  concentrations  of  Cy5.5.  Then 
fluorescence  readings  were  obtained  for  different  source  detector  pairs  and  varying  incident  light 
intensity.  The  acquisition  time  (integration  time)  was  5  sec.  The  concentration  range  examined  was 
[InM  —  800nM]  and  represents  a  wide  span  of  biologically  expected  fluorochrome  accumulations. 
Figure  6  depicts  the  fluorescence  counts  obtained  as  a  function  of  laser  power  for  the  different 
fluorochrome  concentrations  examined.  For  each  light  power  selected,  the  fluorescence  count 
plotted  is  the  actual  fluorescence  measurement  minus  the  background  signal  recorded  by  the  same 
detector  at  OnM  fluorochrome  concentration.  Figures  6a  and  6b  are  the  linear  and  semi-logarithmic 
plots  of  the  fluorescence  counts  obtained  from  the  detector  marked  in  figure  5  with  the  letter  “Y’. 
Similarly,  figures  6c  and  6d  are  the  linear  and  semi-logarithmic  plots  of  the  fluorescence  counts 
obtained  for  the  detector  marked  in  figure  4  with  the  letter  ‘Z’. 

Figure  6  demonstrates  a  hnear  response  for  the  entire  range  of  light  powers  and 
concentrations.  This  finding  indicates  that  no  self-quenching  effects  or  photo-bleaching  occur  with 
the  investigated  indocyanine  dye.  Furthermore  concentrations  of  InM  are  clearly  detectable,  even 
when  using  low  light  power.  This  is  particularly  evident  on  the  semi-logarithmic  plots.  By 
appropriately  selecting  the  light  power  sub-nM  concentrations  could  be  also  detected  for  smaller 
volume  fractions. 

Fluorescence  reconstructions. 

The  feasibility  to  three-dimensionally  reconstruct  and  quantify  fluorochromes  embedded  in 
diffuse  media  with  the  herein  reported  instrument  has  been  demonstrated  recently  In  this  section 
we  investigate  another  aspect  of  the  reconstruction  performance,  i.e.  the  linearity  of  the 
reconstructed  fluorochrome  concentration  as  a  function  of  the  real  fluorochrome  concentration. 
Although  a  linear  relation  between  fluorescence  signals  and  fluorochrome  concentrations  was  shown 
in  Figure  6,  tomographic  techniques  based  on  perturbation  methods  have  been  reported  to  give  non¬ 
linear  responses  and  underestimate  the  perturbation  magnitude  as  this  magnitude  increases.  This  is 
because  perturbation  methods  operate  optimally  on  the  approximation  that  the  perturbation  signal 
detected  is  very  small  compared  to  the  intrinsic  signal  (i.e.  the  signal  at  the  excitation  wavelength). 

To  examine  the  tomographic  performance  we  reconstructed  images  of  the  embedded 
fluorochrome  seen  in  Figiue  5  for  the  range  of  concentrations  examined  in  Figure  6  i.e.  [InM  - 
800nM].  The  laser  power  used  in  all  reconstructions  was  30mW  at  the  fiber  tip  and  the  acquisition 
time  (integration  time)  was  10  seconds  per  source  resulting  in  a  4-minute  acquisition  time  for  a  full 
three-dimensional  reconstruction.  The  reconstruction  segmented  the  field  of  view  in  three  slices  with 
700  voxels  per  slice.  The  voxel  size  was  ~1  x  1  x  8.3  mm^.  The  inversion  technique  used  was  an 
iterative  algebraic  reconstruction  technique  (ART),  which  was  run  until  one  hundred  more  iterative 
steps  yielded  a  change  that  was  less  than  1%.  The  typical  number  of  iterative  steps  was  ~2000. 

A  typical  reconstructed  image  is  seen  in  Figure  7a.  The  image  shown  is  interpolated  to  a  4  x 
4  finer  grid  than  the  original.  Figure  7b  plots  the  reconstructed  chromophore  concentration  as  a 
function  of  the  expected  concentration.  The  reconstructed  value  demonstrates  a  remarkably  linear 
response  for  the  whole  range  of  exantined  concentrations.  It  is  possible  that  perturbation  methods 
apphed  in  the  fluorescence  regime  attain  a  larger  concentration  range  of  hnearity  compared  to 
absorption/scattering  reconstructions.  This  could  be  due  to  the  low  fluorochrome  quantum  yield 
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that  makes  fluotescent  perturbation  signals  much  smaller  than  the  corresponding  absorption  or 
scattering  perturbation  signals. 


DISCUSSION 

We  have  developed  an  imager  for  three-dimensional  reconstructions  of  the  distribution  of 
fluorochromes  in  diffuse  media  and  presented  key  performance  characteristics.  The  detection  fiber 
array  coupling  onto  the  CCD  via  a  lens  system  allows  great  flexibility  in  employing  a  large  number  of 
parallel  channels.  A  set  of  144  fibers  can  be  easily  added  with  the  same  coupling  characteristics  as  the 
36  fibers  investigated  here.  An  arbitrary  large  number  of  fibers  can  be  added  by  constructing  a  larger 
detection  fiber  array.  Consequendy  however  the  field  of  view  acquired  by  the  lens  system  should 
increase  which  will  induce  additional  coupling  losses.  Generally,  a  balance  between  the  detector  array 
si2e  and  the  coupling  efficiency  should  be  achieved  to  attain  an  optimum  design  for  any  given 
application. 

A  particularly  attractive  feature  of  the  fluorescent  tomographer  is  that  sub-nanomolar 
concentrations  of  fluorescent  dyes  contained  into  tumor-sized  volumes  can  be  detected.  The  present 
scanner  was  found  sufficient  to  detect  InM  of  Cy5.5  contained  in  ~100  □  L.  Lower  concentrations  or 
volumes  could  be  detected  by  increasing  light  intensity  and  acquisition  times.  More  sensitive  cameras 
with  lower  noise  floors  (up  to  2  orders  of  magnitude  for  liquid  nitrogen  cooled  CCD  cameras)  could 
also  be  used  to  further  improve  sensitivity.  In  general,  the  sensitivity  levels  achieved  approach  those 
of  nuclear  imaging  techniques  and  are  several  orders  of  magnitude  lower  than  the  sensitivity  for 
detection  of  magnetic  probes  by  Magnetic  Resonance  Imaging  or  x-ray  absorbers  by  Computed 
Tomography. 

We  have  further  demonstrated  that  the  instrument  can  produce  quantified  images  and  that  it 
has  a  linear  response  for  over  two  orders  of  magnitude  of  biologically  relevant  concentrations. 
Generally,  the  linearity  across  the  board  of  parameters  examined  is  an  important  finding  since  no 
correction  schemes  are  requited  and  high  quantification  accuracy  can  be  achieved. 

The  technology  presented  here  is  geared  towards  examining  animals  implanted  with  certain 
diseases  such  as  tumors,  especially  as  pertaining  to  the  bio-distribution  of  targeted  and  activatable 
fluotochrome  reporter  probes  ".  The  specific  geometry  and  dimensions  are  optimized  for  mouse 
imaging.  Accommodation  of  different  size  animals  can  be  achieved  by  using  a  different  size  optical 
bore.  Although  a  large  optical  bore  could  in  principle  be  used  for  smaller  animals  as  well,  optimizing 
the  bore  size  for  the  atumal  size  of  interest  facilitates  minimum  diffusion  and  optimal  photon 
coupling. 

Imaging  of  living  tissues  may  impose  complexities  since  different  background  absorption  or 
scattering  and  spatial  variations  in  the  background  optical  properties  may  result  in  signals  that  are 
different  from  the  ones  obtained  in  our  bench  studies  using  simple  homogeneous  media.  However, 
the  reconstruction  algorithm  employed  is  specifically  designed  for  tissue  measurements.  This  is 
because  the  algorithm  provides  for  a  differential  scheme  where  the  intrinsic  light  normalizes  the 
fluorescence  measurement  to  cancel  out  or  miiumize  the  sensitivity  to  local  and  global  gain  factors 
and  to  heterogeneous  photon  propagation.  An  added  advantage  is  that  this  normalization 
measurement  can  be  obtained  virtually  simultaneously  with  the  fluorescence  measurement  so  that  no 
recordings  before  fluorochrome  administration  are  required.  The  performance  of  this  imager  with 
transgenic  mouse  models,  and  the  cross-examination  of  the  imaging  results  with  other  imaging 
techniques  is  currently  under  investigation. 


FIGURE  CAPTIONS 

Figure  1.  The  animal  tomographer  constmcted. 
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Figure  2.  Optical  bore.  The  sources  and  detectors  are  arranged  in  different  rings  to  allow 
for  three-dimensional  reconstructions. 

Figure  3.  Experimental  signal  to  noise  ratio  of  the  instrument.  The  right  part  of  the 
histogram  depicts  the  signal  to  noise  ratio  achieved  for  different  source  detector  pairs  when  the 
reference  measurement  is  used  for  data  post-processing.  The  right  part  of  the  histogram  plots  the 
signal  to  noise  achieved  when  no  correction  is  applied. 

Figure  4.  a)  Experimental  set-up  that  investigated  the  instrument’s  resolution.  The  two 
objects  shown  are  1.5mm  capillary  tubes,  3  mm  apart,  containing  800nM  of  Cy5.5  dye.  b) 
Reconstruction  achieved  with  the  fluorochromes  inserted  at  4%  intrahpid  (□  s’=5cm  ’).  c) 
Reconstruction  achieved  using  8%  intrahpid  (□  s’=5cm-’). 

Figure  5.  Experimental  set-up  that  used  to  investigate  the  fluorescence  responses  and  the 
reconstruction  hnearity  of  the  technique. 

Figure  6.  Fluorescence  responses  from  the  experimental  set-up  plotted  in  figure  5.  The  top 
line  depicts  in  linear  and  semi-logarithmic  plots  the  responses  detected  by  detector  Y  (c.f  figure  5). 
Bottom  line  depicts  in  linear  and  semi-logarithmic  plots  the  responses  detected  by  detector  Y  (c.f 
figure  5). 

Figure  7.  a)  Fluorescence  reconstruction  of  the  experimental  arrangement  in  figure  5.  The 
figure  shown  is  the  top  shce  of  a  three-shce  reconstmction.  b)  Reconstructed  fluorochrome 
concentrations  as  a  function  of  expected  fluorochrome  concentrations. 
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Abstract 

We  present  a  normalized  Born  expansion  that  facilitates  fluorescence 
reconstructions  in  turbid,  tissue-like  media.  The  algorithm  can  be  particularly  useful  for 
tissue  investigations  of  fluorochrome  distribution  in-vivo  since  it  does  not  require 
absolute  photon-field  measurements  or  measurements  before  contrast  agent 
administration.  This  unique  advantage  can  be  achieved  only  in  fluorescence  mode.  We 
used  this  algorithm  to  three-dimensionally  image  and  quantify  an  indocyanine 
fluorochrome  phantom,  using  a  novel  fluorescence  tomographic  imager  developed  for 
animals. 

Three  dimensional  tomography  and  quantification  of  fluorochrome  bio¬ 
distribution  in  deep  tissue  will  soon  be  necessary  for  imaging  novel  Indocyanine 
derivatives',  targeted  dyes^  and  smart  molecular  beacon-like  agents^  that  enable 
functional  and  molecular  characterization  of  normal  and  diseased  tissues.  Tomographic 
approaches  for  fluorescence  imaging  in  the  NIR  have  been  developed  in  the  past'*’^’®  and 
are  straightforward  expansions  of  algorithms  originally  developed  for  diffuse  optical 
tomography  (DOT)  of  absorption  and  scattering.  In  this  Letter  we  present  an  algorithm 
that  offers  significant  experimental  advantages  because  it  can  work  with  measurements 
taken  only  from  the  target  medium  after  fluorochrome  administration;  no  pre¬ 
administration  measurements  or  phantom  measurements  are  required.  This  could  be 
essential  for  in-vivo  measurements  and  is  a  unique  advantage  of  DOT  in  fluorescent 
mode.  The  algorithm  is  demonstrated  using  an  analytical,  linear  solution  of  the  diffusion 
equation  but  the  principle  of  operation  can  apply  to  non-linear  and  numerical  solutions. 
Experimentally,  we  found  this  algorithm  very  robust  in  imaging  and  quantifying  the 
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three-dimensional  distribution  of  two  fluorochromes  of  different  concentrations 
embedded  into  tissue-like  media  using  a  novel  homemade  animal  imager. 

The  algorithm  developed  is  based  on  a  previously  reported  mainframe  that  uses  a 
Bom-Wpe  approximation  to  predict  the  forward  field  of  photon  propagation  in  diffuse 
media”.  The  photon  wave  detected  in  a  non-bounded  diffuse  medium  at  position  F  due  to 
an  intensity-modulated  source  at  position  F^  is’  Ug(F^,r,(o)  =  0j(FJexp(/AF)/4;cDF, 

where  k  -  +  io))/D)'^^  is  the  photon  wave  propagation  vector,  u  is  the  speed  of 

light  into  the  medium,  o)  is  the  angular  modulation  frequency,  ju^  is  the  absorption 
coefficient,  D  =  u/3/jI  is  the  diffusion  coefficient,  //'  is  the  reduced  scattering 
coefficient  and  (F^ )  is  a  gain  factor  associated  with  light  source  strength,  fiber 

coupling  losses  and  other  attenuation  in  the  optical  system.  If  ©^(F^)  accounts  for  the 
detector  gain  and  the  detection  fiber  losses  and  is  the  detector  quantum  efficiency  at 
wavelength  Xi,  then  the  incident  photon  field  detected  at  position  F^  is: 


where  denotes  the  wave  propagation  vector  for  the  optical  properties  of  the  medium 
of  investigation  at  wavelength  Xi. 

When  fluorescent  objects  exist  in  the  medium,  the  field  f7^(F^,Fj)  detected  at 
position  Pj  at  the  emission  wavelength  can  be  written  as^’*: 


t/^(F„FJ=  pV-0,(FJ.©^  -©.(FJ-C/ofF, 


«(F) 


ia)T{r)  D 


A2 


(2) 


where  n(F)  is  the  product  of  the  fluorochrome  absorption  coefficient  and  fluorescence 
quantum  yield,  r(F)  is  the  fluorescence  lifetime,  k^^  is  the  wave  propagation  vector  at 
wavelength  X.2,  ©y^  is  the  attenuation  of  the  filter  that  is  used  to  collect  the  fluorescent 
field  and  is  the  detector  quantum  efficiency  at  the  emission  wavelength  . 

G(fj  -  r,k^^)  is  the  Green’s  function  solution  to  the  diffusion  equation  and  describes  the 
propagation  of  the  emission  photon  wave  from  the  fluorochrome  to  the  detector.  Eq.2 
treats  fluorochromes  as  two-level  quantum  systems  and  ignores  fluorescence  saturation 
effects  and  high-order  scattering  interactions.  These  approximations  hold  for  the  weak 
fluorescent  concentrations  expected  in  many  biological  systems. 

Finding  solutions  for  Eq.2  requires  determination  of  the  factors  ©i(FJ,©^(F^) 
for  each  source-detector  pair  r,,rj  .To  reach  more  manageable  expressions  for  the 

fluorescent  field  we  can  divide  Eq.2  by  Eq.l  and  rearrange  to  obtain  the  normalized 
fluorescence  Born  field  ,  i.e: 
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Eq.3  is  similar  to  expressions  obtained  for  the  linear  perturbation  DOT  problem 
using  the  Rytov  approximation^’'*’  and  offers  experimental  advantages  because  all  the 
position  dependent  gain  factors  are  cancelled  out.  The  terms  0^  and  QE^^IqE^'  can  be 
determined  experimentally  and  are  virtually  independent  of  position  and  time.  Usually, 
QE^^IqE^'  «!  due  to  the  proximity  of  A,i,  X2. 

Eq.3  assumes  a  measurement  of  the  incident  field  i.e.  a  measurement  through 
a  homogeneous  medium.  This  measurement  is  not  practical  when  probing  biological 
media.  However,  it  has  been  shown"  that  Uinc  can  be  substituted  by  a  measurement  C/^^, 
obtained  through  the  heterogeneous  medium  if  at  the  right  part  of  Eq.  3  is  computed 
for  the  average  optical  properties  of  the  medium  at  (assuming  a  random,  weak  optical 
background  heterogeneity).  It  is  important  to  note  that  can  be  obtained  even  if  the 
fluorochrome  is  present  during  the  measurement.  The  weak  sensitivity  of  the  algorithm  to 
the  selection  of  Ui„c  vs.  U'„^.  is  demonstrated  in  the  experimental  section  of  this  work  and 
is  a  significant  and  unique  advantage  of  the  fluorescence-mode  DOT  compared  to 
differential  absorption/scattering  DOT.  This  is  because  DOT  of  absorption/scattering 
requires  a  measurement  before  fluorochrome  administration'^,  which  may  not  be  always 
available,  due  to  the  long  circulation  times  required  by  some  fluorochromes  for  proper 
target  accumulation. 

For  reconstruction  purposes  Eq.  3  is  discretized  into  a  number  of  volume  elements 
(voxels),  which  yields  a  set  of  coupled  linear  equations.  The  terms  Uo  and  G  can  be 
calculated  numerically  or  analytically  for  the  appropriate  boundary  conditions  and 
background  optical  properties  at  Xi,  X2  respectively.  The  resulting  system  of  equations 
can  then  be  inverted^  for  the  unknown  quantities  n(f)  and  r(F) . 

To  demonstrate  the  algorithmic  performance  we  used  a  novel  animal  optical 
tomographer  with  24  source  and  36  detector  fibers  arranged  in  five  equidistant  rings 
around  a  cylinder  as  shown  in  Figure  1.  The  cylinder  was  filled  with  0.5%  Intralipid 
(Fresenius  Kabi  Clayotn  L.P.,  NC)  and  6pM  IC-Green  (Akom  Inc.  IL)  yielding  a  liquid 
with  //'=7  cm''  and  =0.05  cm''  at  675nm.  Two  glass  tubes  (4mm  outer  -  3.5  mm 
inner  diameter)  were  inserted  with  centers  12  mm  apart,  at  different  depths  as  shown  in 
Figure  1.  The  first  tube  that  traversed  the  field  of  view  FOV  (i.e.  the  volume 
reconstructed)  contained  500nM  of  Cy5.5  (Amersham  pharmacia  biotech  NJ)  and  0.5% 
Intralipid  and  the  second  tube  that  was  inserted  only  by  8mm  in  the  FOV  contained 
750nM  of  Cy5.5.  Cy  5.5  concentrations  were  verified  with  a  calibrated  U3000  Hitachi 
spectrophotometer.  The  light  source  was  a  675nm,  2m W  continuous  wave  laser  diode 
(Thorlabs  NJ)  with  2%  amplitude  noise.  Light  was  switched  sequentially  to  the  24  source 
positions  using  an  optical  switch  (Dicon  FiberOptics  CA).  The  light  was  directed  onto  the 
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optical  bore  using  a  200/230|j,m  multimode  optical  fibers  (Fiber  Instrument  Sales  NY). 
For  light  detection  3mm  diameter  optical  fiber  bundles  were  used  (Dollan-Jenner 
Industries  MA)  which  were  imaged  on  a  CCD  camera  (Roper  Scientific  NJ)  using  a 
Micro-Nikkor  f/2.8D  AF  lens  (Nicon  JP).  For  fluorescence  detection  a  three-cavity  band¬ 
pass  filter  was  used  with  center  frequency  at  710nm  and  lOnm  full- width  at  half 
maximum  (Andover  Corporation  NH). 

Inversion  of  Eq.  3  after  discretization  was  performed  using  the  algebraic 
reconstruction  technique®,  which  was  run  until  the  result  showed  changes  of  less  than 
0.01%  per  10  iterations  (approximately  3000  iterations  run  in  ~  5  minutes  on  a  IGHz 
Intel  Pentium  III  personal  computer).  Since  the  problem  is  solved  at  ccf=0,  only  n(r)  was 
reconstructed.  The  functions  Uo  and  G  were  calculated  for  the  cylindrical  boundary 
using  a  finite-difference  solution  of  the  homogenous  diffusion  equation  for  the 
background  medium.  A  zero-boundary  condition'^  was  assumed  at  1/ //'  outside  the  actual 

boundary.  Three  planes  were  simultaneously  reconstructed  along  the  z  plane  with  voxel 
size  /i=0.833x0.12x0.12  cm^  (total  of  1047  voxels). 

Figure  2  depicts  the  reconstructed  planes  using  Ui„c  in  Eq.3,  which  is  the  incident 
photon  field  measured  through  the  homogeneous  background  solution  of  IC-Green  and 
Intralipid  with  no  fluorochrome  objects  immersed.  Each  plane  depicted,  is  interpolated  to 
a  60x60  grid  (h  -  0.04x0.04  cm^).  Two  cylinders  appear  in  the  top  plane  and  one  cylinder 
in  the  middle  and  lower  plane  as  expected.  Quantification  was  within  10%  of  the 
expected  Cy5.5  concentration.  The  size  of  the  cylinder  was  correctly  resolved  at  the 
middle  plane  (within  ±  half  voxel  size)  but  was  overestimated  for  the  top  and  bottom 
planes  by  ~30%.  This  is  typical  in  such  reconstructions  because  the  top  and  bottom 
planes  have  virtually  half  the  number  of  light  path-length  projections  passing  through 
them  than  the  planes  in  the  middle  and  hence  there  is  less  information  available  for  the 
medium  at  the  top  and  bottom.  To  overcome  this  effect  it  would  be  desirable  to  increase 
the  number  of  source-detector  rings  along  z  so  that  the  field  of  view  is  fully  contained 
within  the  source-detector  rings. 

Figure  3  depicts  the  reconstructions  obtained  when  using  in  Eq.3,  which  is 
the  photon  field  measured  at  the  excitation  wavelength  with  the  fluorochrome  objects 
immersed.  Although  no  short-wave-pass  filter  was  used,  the  influence  of  emission 
photons  is  negligible  because  U'„c»Uj, .  The  reconstructed  images  appear  essentially 

identical  to  the  ones  in  Figure  2.  There  were  small  (<3%)  differences  in  the  quantification 
of  objects.  The  shape  of  the  objects  reconstructed  also  demonstrates  minor  differences 
but  the  center  of  the  objects  did  not  change.  Tissue-like  background  heterogeneity  could 
have  a  stronger  effect  in  the  accuracy  of  the  result.  However  this  effect  should  be  no 
different  than  the  one  seen  in  clinical  applications  of  DOT  for  imaging  absorption 
perturbations’^,  since  absorption  DOT  uses  analogous  assumptions". 

In  conclusion,  we  have  demonstrated  three-dimensional  tomographic 
reconstructions  of  different  fluorochrome  concentration  in  diffuse,  tissue-like  phantoms 
using  a  novel  animal  imager.  These  results  were  obtained  using  a  reconstruction 
algorithm  that  is  appropriate  for  clinical  measurements  since  measurement  calibration  can 
be  performed  virtually  simultaneously  with  the  fluorescent  measurement,  hence 
maintaining  high  positional  accuracy,  under  the  same  physiological/biochemical 
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conditions  and  maximize  patient  comfort.  The  algorithm  and  described  tomographic 
imaging  system  should  provide  a  highly  useful  tool  to  measure  a  variety  of  functional  and 
molecular  abnormalities  in  tissues  in-vivo. 


Figure  1.  a)  Optical  bore,  made  of  Delrin®  acetal  resin.  The  fibers  pass  through  the  Delrin 
material  and  come  face  with  the  inner  walls.  The  tubes  immersed  into  the  intralipid 
solution  are  shown  in  dotted  lines,  b)  Side  view  of  (a). 
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Figure  2.  Fluorescence  reconstruction  of  the  experimental  set-up  of  Figure  1  using  three 
layers  along  z.  The  fluorescence  measurements  are  normalized  using  Umc  i.e.  the 
incident  field  measured  when  the  two  absorbing/fluorescent  objects  were  not  immersed. 
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Figure  3.  Fluorescence  reconstruction  of  the  experimental  set-up  of  Figure  lb  using  three 
layers  along  z.  The  fluorescence  measurements  are  normalized  using  U’mc  i.e.  the 
incident  field  measured  with  the  two  absorbing/fluorescent  objects  immersed. 
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Abstract 

Systematic  efforts  are  under  way  to  develop  novel  technologies  that  would  allow 
molecular  sensing  in  intact  organisms  in  vivo.  Using  inversion  techniques  that  account  for 
the  diffuse  nature  of  photon  propagation  in  tissue  and  near  infrared  fluorescent  molecular 
beacons  we  were  able  to  obtain  three-dimensional  in-vivo  images  of  cathepsin  B  expression 
of  orthopic  gliomas.  We  demonstrate  that  activatable  fluorochromes  can  be  detected  with 
high  positional  accuracy  in  deep  tissues,  that  molecular  specificities  of  different  beacons 
towards  enzymes  can  be  resolved,  and  that  tomography  of  beacon  activation  is  linearly 
related  to  enzyme  concentration.  The  tomographic  imaging  method  offers  a  range  of  new 
capabilities  for  studying  biological  function  using  fluorescent  chemical  sensors,  for 
identifying  molecular  expression  patterns  via  multispectral  imaging  and  for  continuously 
moititoring  drug  therapies. 

Keywords:  molecular  sensing,  murine  models,  tomography,  imaging,  near  infrared 
fluorescence 
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Main  text 


Compounds  exhibiting  fluorescent  properties  in  the  visible  and  near-infrared  Ught 
range  play  an  important  role  in  many  biotechnology  applications  including  gene  expression 
profiling,  determining  protein  function,  elucidating  cellular  pathways,  in  sensing  small 
molecule  protein  interaction,  or  in  measuring  physiologic  solutes  Inducible 
photoproteins,  sensing  protein  chimera,  organic  fluorochromes,  quantum  dots,  resonating 
nanoparticles  and  lanthanide  complexes  have  thus  been  used  as  reporters  for  probing 
fluorescence,  fluorescence  energy  transfer  (FRET),  spectral  analysis  or  fluorescence  lifetime 
More  recently,  near  infrared  (700-900  nm)  targeted  fluorochromes,  sensing 
fluorogenic  substrates  and  red-shifted  photoproteins  with  higher  tissue  penetration  have 
been  developed  to  measure  molecular  events  in  vivo  in  small  animals 

Fluorescence  generated  in  in-vitro  assays  can  be  easily  quantified  using  fluorometers 
or  charge  coupled  devices  (CCD).  Similarly,  fluorescence  of  superficial  stmctures  has  been 
imaged  in  vivo  using  confocal  or  multiphoton  microscopy  Quantitation  of  fluorescence 
by  imaging  and  sensing  fluorochrome  reporters  in  deep  stmctures  (>  1  mm)  however,  has 
been  more  elusive.  Here  we  describe  the  in-vivo  implementation  of  fluorescence  mediated 
tomography  (FMT)  to  image  fluorescence  in  deeper  tissues  in  murine  models  and  show  how 
highly  sensitive  measurements  with  molecular  specificity  can  be  obtained.  FMT  is  a  novel, 
quantitative  3D  imaging  technique  capable  of  sensing  pico-  to  femtomole  quantities  of 
fluorochromes  in  tissues  at  macroscopic  scale,  ie.  in  whole  animals  with  mm  resolution.  The 
technique  shares  tomographic  principles  with  diffuse  optical  tomography  but 
simultaneously  uses  absorption  and  fluorescent  measurements  for  accurate  three- 
dimensional  reconstmction  of  fluorochrome  concentration 

A  tomographic  imaging  system,  shown  in  Fig.  1,  was  constmcted.  Laser  light  at  670 
nm  (a),  operating  in  40  mW  continuous  mode,  was  used  to  successively  excite  24  point 
sources  in  two  12-fiber  layers  within  the  cylindrical  imaging  chamber  (e)  (2.5  cm  diameter) 
by  using  an  optical  switch  (d)  (Dicon  Fiberoptics  CA).  For  each  excitation  point,  the  output 
of  three  layers,  each  of  12  detection  fiber  bundles  was  captured  onto  a  6x6  array  (h)  (c.f. 
inset  of  Fig.l)  and  imaged  using  a  high  quantum  yield  CCD  camera  (j)  (Roper  Scientific  NJ). 
Measurements  were  acquired  for  excitation  light  (absorption  and  scattering  signal),  and 
fluorescence  using  a  band-pass  interference  filter  (i)  (Andover  Corporation  NH)  yielding  two 
vectors  of  864  data  points  each.  Additional  measurements  were  obtained  in  a  reference 
channel  (c)  through  a  beam-splitter  (b)  (OZ  Optics  Ontario  Canada)  to  account  for  temporal 
laser  fluctuations.  For  image  reconstmction  and  data  processing,  the  illuminated  imaging 
volume  was  segmented  into  5  layers  containing  each  491  voxels,  the  voxel  size  being  1  x  1  x 
5  mm^  For  data  reconstmction,  each  voxel  was  assigned  an  initially  unknown  fluorochrome 
concentration.  Then  a  forward  problem  was  constmcted  that  predicted  a  synthetic 
measurement i.e.  the  ratio  of  photons  detected  at  the  emission  wavelength  after  correction 
for  bleed-through  signals  over  the  photons  detected  at  the  excitation  wavelength  for  an 
assumed  known  fluorescence  distribution  in  the  medium  (c.f.  inset  of  Fig.l)  .  The  forward 
model  was  inverted  using  algebraic  reconstmction  techniques  with  positive  restriction  for 
the  experimentally  obtained  measurement  vector.  Reconstmcted  images  display  the 
fluorochrome  concentration  in  nanomolar  units,  calibrated  against  free  and  quenched  near 
infrared  indocyanine  compounds. 
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To  determine  the  sensitivity  by  which  measurements  could  be  carried  out  in  turbid 
media,  indocyanine  fluorochrome  samples  (100  pL)  were  embedded  in  highly  scattering 
media  simulating  optical  properties  of  parenchymal  tissue.  With  the  current  acquisition  set¬ 
up,  the  detection  threshold  of  fluorochrome  was  approximately  two  hundred  femtomoles  in 
the  embedded  target  volume  deep  in  the  phantom  The  detection  limit  is  a  function  of 
multiple  variables  including  the  number  of  source/ detector  pairs,  imaging  time,  laser  output, 
CCD  characteristics  and  the  optical  properties  and  dimensions  of  the  medium  of 
investigation.  With  state  of  the  art  technology  it  is  estimated  that  fluorescence  detection 
could  be  further  decreased  by  at  least  two  orders  of  magnimde  using  liquid  nitrogen  cooled 
CCD  cameras,  higher  power  lasers,  more  red-shifted  fluorochromes  and  by  improving  the 
coupling  efficiency  of  photons  onto  the  CCD  chip  In  a  separate  experiment  the  spatial 
resolution  of  the  constmcted  tomographer  was  also  determined  by  the  ability  to  separate  two 
1  mm  cylinders  embedded  in  fluid  simulating  optical  properties  of  parenchymal  tissue.  The 
spatial  resolution  of  the  FMT  system,  ranged  from  1  mm  near  the  surface  to  approximately 
2.5  mm  in  the  center  of  the  volume  Finally  the  quantification  error  has  been  shown  to  be 
within  10%  of  the  expected  fluorochrome  concentration  for  the  experimental  set-up 
employed  Again,  increasing  the  source/detector  pairs  and  power  settings  could  be  further 
used  to  improve  the  above  resolutions  and  quantifications  by  improving  the  spatial 
information  and  signal  to  noise  ratio  detected. 

The  ability  to  detect  fluorochromes  in  vivo  would  be  of  considerable  biomedical 
interest,  as  generic  fluorochromes  and  those  targeted  to  receptors  have  recently  been 
described  An  inherently  higher  signal-to-noise  ratio  can  be  achieved  by  using  activatable 
fluorochromes  (“molecular  beacons”),  which  become  fluorescent  after  hybdridization  or 
enzymatic  cleavage  We  therefore  tested  whether  the  tomographic  imaging  system  could 
resolve  the  enzyme  activity  of  cathepsin-B  in  turbid  media  using  a  previously  developed 
biocompatible  NIR  cathepsin  B  sensitive  molecular  beacon  The  experimental  setup, 
shown  in  Fig  2a,  was  such  that  four  capillaries  contained  200  pL  of  a  1  pM  cathepsin-B 
sensitive  imaging  probe  to  which  was  added  either  0,  5,  25  or  50  pg  of  purified  cathepsin  B 
(Calbiochem,  La  Jolla,  CA).  The  capillaries  were  enclosed  in  turbid  media  as  shown  in  Fig.  2a 
and  imaged  over  time  with  each  acquisition  taking  2  minutes.  The  surrounding  turbid  media 
was  constructed  out  of  polyester  resin  (Creative  Wholesale  Stockbridge  GA)  to  which  Ti02 
particles  (Sigma  St.  Louis  MI)  and  India  Ink  (KOH-I-NOOR  Leeds,  MA)  had  been  added  to 
simulate  tissue  scattering  and  absorption.  The  reconstructed  values  shown  in  Fig.  2b 
demonstrate  a  linear  relation  between  enzyme  concentration  and  fluorochrome  activation 
covering  a  range  of  physiologically  encountered  enzyme  levels,  obtained  at  4  and  24  hours 
after  enzyme  administration.  This  result  indicates  that  fluorescence  reconstmctions  could  be 
used  for  quantitative  insights  into  enzymatic  function. 

In  another  set  of  experiments,  we  examined  whether  the  tomographic  imaging 
system  could  resolve  specific  enzymatic  activity  when  tested  with  the  above  described 
reporter  probe.  Cathepsin  B,  D  and  H  were  added  in  different  capillaries  of  the  same 
experimental  setup  shown  in  Fig.  2a  at  25  pg/capiUary  diluted  in  buffer  AU  capillaries 
contained  the  cathepsin  B  sensitive  NIR  fluorescent  probe  (1  pM).  The  fourth  capillary  was 
in  the  absence  of  an  enzyme  for  control  purposes.  Fig.  2c  shows  the  reconstructed  middle 
slice  24h  after  substrate/enzyme  interaction.  Activation  of  the  cathepsin  B  sensitive  probe 
by  cathepsin  B  is  clearly  evident,  whereas  it  was  lacking  for  the  other  cathepsins.  The  kinetics 
of  probe  activation  due  to  cathepsin  B  versus  cathepsin  D  is  shown  in  Fig  2d.  Cathepsin  H 
demonstrated  a  very  similar  activation  pattern  with  cathepsin  D  while  the  fourth  tube 
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demonstrated  minimal  background  fluorescence.  Measurements  performed  under  identical 
conditions  in  non-turbid  media  using  a  fluorometer  (F4500,  Hitachi)  are  plotted  on  Fig.  2d 
as  triangles  and  verify  the  activation  rate  calculated  by  FMT.  Summarily,  the  above  results 
show  that  1)  the  technique  is  specific  in  differentiating  fluorescence  activation  imparted  by 
different  enzymes,  that  2)  there  is  a  predictable  relation  between  enzyme  activity  and 
generated  fluorescence  under  the  tested  conditions  and  that  3)  measurements  obtained  by 
FMT  imaging  from  within  turbid  media  closely  reflect  those  obtained  using  an  analytical 
fluorometer  under  non-turbid  conditions.  It  appears,  that  the  technique  should  thus  be  able 
to  measure  a  number  of  different  sensing  molecular  probes  currently  under  development. 

While  the  above  experiments  were  performed  in  vitro  using  phantoms  and  tissue-like 
turbid  media,  the  main  challenge  was  to  determine  whether  similar  experiments  could  be 
performed  in  murine  models.  As  an  extension  to  the  above  validation  studies,  we  chose  to 
image  cathepsin  B  activity  in  9L  gliosarcomas  stereotactically  implanted  into  unilateral  brain 
hemispheres  of  nude  mice,  as  cathespin  B  activity  had  been  implicated  in  glioma  invasion 
26,29,30  ^  animals  were  subjected  to  correlative  MR  imaging  to  determine  the  presence  and 
location  of  tumors  prior  to  the  FMT  imaging  studies.  Fig.  3  summarizes  a  representative 
imaging  experiment.  Fig.3a-b  depict  the  gadolinium-enhanced  tumor  (enhancement  is  shown 
in  a  green  color  map  superimposed  onto  a  T1  weighted  image)  on  axial  (a)  and  sagital  (b) 
slices.  Fig.  3c,  e  and  f  depict  the  three  consecutive  FMT  slices  obtained  from  top  to  bottom 
of  the  volume  of  interest.  The  location  and  volume  covered  by  the  three  slices,  is  indicated 
on  Fig.  3b  by  thin  white  horizontal  lines.  Fig.  3c  shows  marked  local  probe  activation  relative 
to  adjacent  slices,  congruent  with  the  location  of  the  tumor  identified  on  the  MR  images. 
Fig.  3d  shows  a  superposition  of  the  MR  axial  slice  passing  through  the  mmor  (Fig.3a)  onto 
the  corresponding  FMT  shce  (Fig.  3c)  after  appropriately  translating  the  MR  image  to  the 
actual  dimensions  of  the  FMT  image.  For  coregistration  purposes  we  used  special  water- 
containing  fiducials  and  body  marks  that  facilitated  matching  of  the  MR  and  FMT 
orientation  and  animal  positioning.  The  in-vivo  imaging  data  correlated  well  with  surface 
weighted  reflectance  imaging  of  the  excised  brain.  Fig.  3g  depicts  the  axial  brain  section 
through  the  9L  tumor  examined  with  white  light  using  a  CCD  camera  mounted  onto  a 
dissecting  microscope.  Fig.  3h  shows  the  same  brain  section  imaged  using  a  previously 
developed  reflectance  imaging  system  at  the  excitation  wavelength  (675  nm)  and  Fig.  3i  is 
the  fluorescence  image  obtained  at  the  emission  wavelength  using  appropriate  3-cavity  cut¬ 
off  filters  (Andover  Corporation  NH)  that  demonstrates  a  marked  fluorescent  probe 
activation,  congment  with  the  tumor  position  identified  by  gadolinium-enhanced  MRI  and 
FMT.  On  average,  the  concentration  of  activated  fluorochrome  per  tumor,  as  calculated  by 
FMT  was  about  110±30  nM.  The  above  results,  confirm  that  cathepsin  B  can  be  used  as  an 
imaging  marker  since  the  protease  is  produced  in  considerable  amounts  by  tumor  cells  and 
recruited  host  cells  Cathepsin  B  expression  in  the  tumors  was  further  confirmed  by 
immunohistochemistry.  Western  blotting  and  RT-PCR  (data  not  shown). 

To  study  whether  FMT  imaging  could  reveal  different  levels  of  cathepsin  B  enzyme 
expression  we  also  performed  correlative  imaging  and  analysis  in  another  tumor  model 
(HT1080)  known  to  have  higher  levels  of  a  number  of  different  proteases  Using  the 
same  imaging  set-up  described  above  we  determined  a  5  fold  higher  cathepsin  B  activity  in 
size  matched  HT1080  tumors  (a  human  fibrosarcoma)  compared  to  9L  tumors.  These 
results  were  confirmed  by  RT-PCR  and  western  blotting,  the  latter  revealing  4  fold  higher 
cathepsin  B  protein  levels  in  HT1080  tumors  compared  to  9L  tumors.  Similar  results  were 
also  observed  in  other  models  in  which  HT1080  tumors  had  been  implanted  into  the  lung  or 
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retroperitoneum  or  when  transgenic  animals  developing  breast  tumors  were  used  (data  not 
shown).  FMT  imaging  thus  offers  the  potential  for  investigating  enzyme  specific  events  in 
live  animals,  opening  unique  opportunities  for  sensing  a  wider  array  of  enzymes.  While  the 
current  studies  only  used  cathepsin  B  as  a  marker  for  consistency,  many  other  proteases 
could  potentially  be  revealed  in  vivo.  In  more  recent  experiments  biocompatible  imaging 
probes  have  been  synthesized  for  cathepsin  D  cathepsin  K,  MMP-2  among  others.  The 
number  of  known  proteases  is  increasing  at  a  tremendous  rate  as  a  consequence  of  genome 
sequencing  projects.  Through  sequence  homology  to  known  proteases,  and  other  high 
throughput  tools  it  should  be  possible  to  find  unique  peptide  substrates  with  high 
specificity  for  given  enzymes.  Furthermore,  larger  numbers  of  affinity  molecules,  derived 
from  high  throughput  screens  could  be  ta^ed  with  fluorochromes  and  be  rapidly  screened 
in  mouse  models.  Such  “in  vivo  assays”  could  become  powerful  adjuncts  to  dmg 
development  as  they  take  into  account  physiology  and  disease-host  interactions.  With  image 
acquisition  times  of  only  a  few  minutes  (typically  two  to  three  minutes  in  the  animals 
studied),  FMT  imaging  has  considerable  advantages  over  existing  technologies  in  providing 
high  throughput  functional  information  in  intact  living  animals. 

How  can  the  technology  be  adapted  clinically?  NIR  light  has  been  shown  to 
penetrate  tissues  several  centimeters  and  detection  is  thus  not  necessarily  limited  to  small 
volumes.  Indeed,  recent  measurements  and  modeling  show  that  near  infrared  fluorescent 
signals  from  tumor-like  stmctures  can  propagate  for  more  than  15  centimeters  in  breast  or 
lung  tissue  and  more  than  5  cm  in  the  adult  brain.  Again,  these  estimates  may  well  be 
surpassed  by  newer  detection  technologies  and/or  more  efficient  beacons.  While  we  chose 
to  construct  a  cyUndrical  imaging  chamber  for  mice,  human  imaging  systems  could  well  have 
different  geometries  for  more  efficient  light  collection.  Planar  imaging  geometries,  handheld 
scanners  or  even  endoscopic  sensors  are  all  technologically  feasible  providing  optimized 
approaches  for  imaging  of  specific  organs  or  diseases.  An  important  aspect  to  potential 
clinical  tomographic  imaging  is  the  fact  that  FMT  is  inherently  quantitative  and  can  be 
acquired  over  time  for  monitoring  purposes.  Additional  advantages  include  the  fact  that  no 
radiation  is  required,  that  beacons  and  fluorochromes  are  usually  stable  and  do  not  decay  like 
isotopes  and  that  the  technology  is  relatively  inexpensive  compared  to  other  tomographic 
imaging  systems.  It  is  also  conceivable  to  perform  multi-wavelength  imaging  to  obtain 
information  from  multiple  targets  or  to  validate  measurements  similar  as  in  comparative 
hybridization  experiments.  The  former  is  of  particular  clinical  relevance  as  “expression 
patterns”  rather  than  single  molecules  may  form  the  basis  for  more  sophisticated  future 
imaging  techniques.  Finally,  FMT  imaging  should  be  useful  for  monitoring  drug  therapies  at 
the  molecular  levels  ^  over  time  and  in  tailoring  treatments  in  individual  patients.  FMT 
imaging  may  also  be  a  useful  adjunct  to  x-ray  computed  tomography  and  magnetic 
resonance  imaging.  It  is  thus  hoped  that  FMT  imaging  will  find  applications  in  unraveling 
complex  biological  pathways  in  murine  models  and  in  applying  this  knowledge  to  current 
problem  areas  in  clinical  diagnosis  and  drug  testing. 
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FIGURES 


Figure  1.  Schematic  of  the  FMT  imager  used  for  the  experimental  measurements.  The  insert 
surrounded  by  the  dotted  line  describes  the  typical  measurements  and  formation  of  the 
composite  measurement,  necessary  to  produce  quantified  images  (for  details  see  text). 
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Figure  2.  a)  Photograph  of  the  turbid  resin  tube  (shown  inserted  into  the  optical  chamber), 
used  for  two  in-vitro  experiments  summarized  here.  The  four  capillaries  shown  were  filled 
with  IgM  of  a  cathepsin-B  sensitive  NIR  activatable  probe  for  both  experiments,  b) 
(Experiment  1);  Probe  activation  as  a  function  of  cathepsin  B  concentration  obtained  with 
FMT  at  4  and  24  hours  after  enzyme  addition,  c)  (Experiment  2):  FMT  reconstruction 
obtained  from  the  middle  slice  of  the  resin  tube  at  24h  after  25  pg  of  cathepsins  (B,  D,  H) 
were  added  in  different  capillaries  as  marked  on  the  image.  The  fourth  capillary  did  not 
contain  any  enzyme,  d)  (Experiment  2  cont.’):  Activation  kinetics  obtained  with  FMT  from 
the  capillaries  containing  cathepsin  B  and  D.  The  orange  triangles  indicate  measurements 
obtained  in  parallel  with  a  fluorometer  in  the  absence  of  the  turbid  medium. 
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Figure  3.  In-vivo  FMT  of  cathepsin  B  expression  levels  in  9L  gUosarcomas  stereotactically 
implanted  into  unilateral  brain  hemispheres  of  nude  mice,  a),  b)  Axial  and  sagital  MR  slices 
of  an  animal  implanted  with  a  tumor,  which  is  shown  in  green  after  gadolinium 
enhancement,  c),  e),  Consecutive  FMT  slices  obtained  from  top  to  bottom  from  the 
volume  of  interest  shown  on  (b)  by  thin  white  horizontal  lines,  d)  Superposition  of  the  MR 
axial  slice  passing  through  the  tumor  (a)  onto  the  corresponding  FMT  slice  (c)  after 
appropriately  translating  the  MR  image  to  the  actual  dimensions  of  the  FMT  image,  g),  h) 
Axial  brain  section  through  the  9L  tumor  imaged  with  white  light  and  with  monochromatic 
light  at  the  excitation  wavelength  (675  nm)  respectively,  and  i)  fluorescence  image  of  the 
same  axial  brain  section  demonstrating  a  marked  fluorescent  probe  activation,  congruent 
with  the  tumor  position  identified  by  gadolinium-enhanced  MRI  and  FMT. 
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1 9.  The  exact  form  of  the  synthetic  measurement  C/y,  (^ ,  ^  )  detected  at  position  due  to  a  source  at 


position  ^ ,  and  its  relation  to  the  fluorescent  distribution  in  the  tissue  of  investigation  can  be 
written  as 


t/(F,,FJ  =  5o 
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where  Uj, (7^,rj) are  measurements  at  the  emission  and  excitation  wavelength 
respectively,  i!7y  >  O  )  ~  ® /  ‘  (^  >  O  )  '*  bleed  through  signal,  ©y  is  the  band-pass 
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reports  written  for  this  Command.  Request  the  limited 
distribution  statement  for  the  enclosed  accession  numbers  be 
changed  to  "Approved  for  public  release;  distribution  unlimited. 
These  reports  should  be  released  to  the  National  Technical 
Information  Service. 

2.  Point  of  contact  for  this  request  is  Ms.  Kristin  Morrow  at 
DSN  343-7327  or  by  e-mail  at  Kristin.Morrow@det . amedd. army.mil . 

FOR  THE  COMMANDER: 


Deputy  Chief  of  Staff  for 
Information  Management 
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